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Abstract 

In this paper we analyze selected evolution models of A—level open quantum systems in order 
to find the minimal number of observables (Hermitian operators) such that their expectation 
values at some time instants determine the accurate representation of the quantum system. The 
assumption that lies at the foundation of this approach to quantum tomography claims that time 
evolution of an open quantum system can be expressed by the Kossakowski - Linclblad equation 
of the form p = Lp, which is the most general type of Markovian and time-homogeneous master 
equation which preserves trace and positivity. We consider the cases when the generator of 
evolution can be presented by means of two or more double commutators. Determining the 
minimal number of observables required for quantum tomography can be the first step towards 
optimal tomography models for A—level quantum systems. 


1 Introduction 

In this paper we shall denote the Hilbert space by hi and we shall assume that dim'H = N < oo. 
B(7i) shall refer the complex vector space of all bounded linear operators in hi. Naturally, B{hi ) is 
isomorphic with the space of N-dimensional complex matrices. The latter shall be represented by 
Mjv(C). Finally, B*(hi) shall refer to the real Banach space of self-adjoint (hermitian) operators on 
hi. In Physics the elements of B*(hi) are called observables due to the fact that an element from 
this space is assigned to every measurable quantity . 

The term quantum tomography is used in reference to all methods and approaches which aim to 
reconstruct the accurate representation of a quantum system by conducting a series a measurements. 
Among a large number of approaches to quantum tomography one can especially mention the so- 
called static model of tomography, which requires N 2 — 1 measurements of distinct observables 
taken at time instant t = 0 (see more in [1, 2, 3]). Evidently, as the number of observables increases 
with the square of N this approach seems rather impracticable. Therefore, there is a need to devise 
other, more efficient, models - for example in [4] one can read about reconstruction of the density 
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matrix by simple projectors. A paper published in 2011 introduced a new approach to quantum 
tomography which is based upon weak measurement. In that paper it has been proved that the wave 
function of a pure state can be measured directly [5], in a contrast to a common belief. Subsequent 
papers revealed that this approach can be generalized also for mixed state identification [6]. 

In this paper we follow yet another approach to quantum tomography - the stroboscopic tomog¬ 
raphy (or stroboscopic observability) which initiated in 1983 in the article [7]. This approach was 
developed in further papers, such as [8, 9]. One can also refer to a well-written review article [10], 
which contains all fundamental results. The main advantage of this method is that it enables to 
determine the optimal criteria for observabiliy of an open quantum system. Thus it seems to have 
considerable potential applications in experiments. 

The assumption that lies at the very foundation of the stroboscopic tomography claims that the 
evolution of an open quantum system can be expressed by a master equation of the form 


P = T p , (1) 

where the operator L shall be referred to as the generator of evolution and its most general structure 
in diagonalized form can be expressed as [ 11 , 12 ] 

TV 2 —1 

Lp = -i[H, p] + -Y J H {[V iP , Vfl + [Vi,pV*]) , (2) 

2=1 

where H € B^iH) and 7 \ > 0. Operators V) G B(TL ) are called Lindblad operators. 

One can easily notice that if V) = V* for i = 1,... ,N 2 — 1, one is able to rewrite the general 
form of the generator of evolution as 

tv 2 -i 

L p = -i[H,p\-~ £ 7 i\Vi,Wi,p]], (3) 

2=1 

where the dissipative part of the generator is given by means of double commutators. Generators 
of the form given by (3) have been the subject of many analysis in quantum Physics. Recent 
results introduced in [13] show that the double commutator form of the generator is convenient for 
analyzing qubit decoherence in the framework of geometric quantum mechanics. Thus the main 
motivation for the current paper comes from the observation that double commutators appear in 
the structure of the generator of evolution under the assumption that the Lindblad operators are 
Hermitian along with the fact that the generator of structure given by (3) is considered in other 
aspects of quantum Physics. 

In order to determine the initial density matrix p( 0) (and on the basis of the knowledge about 
the evolution - the complete trajectory of the state) one needs to assume the availability of a set of 
identically prepared quantum systems, each with evolution given by the master equation with the 
generator L. Furthermore, one has to bear in mind that each system can be measured only once, 
because any measurement, in general, changes the state. 

In the stroboscopic tomography one assumes to have a set of r distinct hermitian operators 
denoted by {Qi, ... ,Q r } ( Q* = Qi ), where r < N 2 — 1 . Each operator can be measured at one 
or more time instants from the set {ti,... ,tp}- The measurement results, denoted by rrii(tj), are 
represented by the formula rrii(tj ) = Tr(Qip(tj)). From many possible research problems connected 
with the stroboscopic tomography, in this paper we shall discuss the aspect of the minimal number 
of observables required for quantum tomography, which is considered the most important. One can 
recall the theorem [9]. 
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Theorem 1.1. For a quantum system with the time evolution given by the Kossakowski-Lindblad 
master equation of the form (1), there exists a number (denoted byp) which determines the minimal 
number of observables required to reconstruct the density matrix. It can be computed from the equality 

r] := max {dim.KerfL — AI)}. (4) 

ASo-QL) 

The number p shall be called the index of cyclicity of the quantum system. 

If one wants to employ a quantum tomography model in an experiment, then the index of 
cyclicity p tells them how many distinct experimental setups one needs to prepare to be able to 
reconstruct the density matrix. In general, it is more efficient to repeat the same kind of measure¬ 
ment a few times then to prepare a few different kinds of measurements. Therefore, we believe that 
the stroboscopic tomography possesses significant potential for future applications as it focuses on 
determining the optimal criteria for quantum tomography, i.e. the minimal number of distinct ob¬ 
servables. One can observe that the value of p depends on the geometric properties of the generator 
of evolution, but the meaning of p is purely physical. 

The well-known results concerning the index of cyclicity for N— level quantum systems relate to 
the time evolution given by a von Neumann equation (see the result in [7]) and Gaussian semigroups, 
i.e. the generator of evolution given by a double commutator L p = — ^[H, [H, p\] where H £ B*(TT) 
(see more in [9]). In this article we consider more general evolution models when the generator is 
given by two or more double commutators. From mathematics we can adopt methods that allow 
to analyze such operators with respect to the eigenvalues and minimal polynomial (see for example 

[14])- 

This paper shows that under certain assumptions one can successfully calculate the index of 
cyclicity for an IV—level open quantum system for a variety of evolution models. This work can be 
the first step towards complete and optimal tomography models for N— level quantum systems. 

2 Generator of evolution with two Lindblad operators 

In this section we shall analyze an open quantum system with the generator of evolution that 
consists of two Lindblad operators, which shall be denoted by G and F (of course, F,G € B{T~i)). 
Thus the evolution of the quantum system in question is given by 

P = \ ([Fp, F*) + [F, P F *] + [Gp, G*] + [G, P G *]). (5) 

In order to be able to discuss the problem of observability of a system with such evolution we 
have to assume a kind of relation between the Lindblad operators F and G. In this part we propose 
to consider the case when F = G*. With this assumption it is possible to rewrite equation (5) in 
the form 

P = \ ([Fp , F*} + [F, pF*} + [F*p, F] + [. F*,pF ]). (6) 

One can easily observe that the generator of evolution can be expressed by means of two double 
commutators 

Lp = “([ F*,[F,p]} + [F,[F*,p}}). (7) 

If one wants to obtain a specific result concerning the index of cyclicity, it is necessary to make 
more assumptions about the operator F and its eigenvalues. First, let us denote by a(F ) the 
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spectrum of the operator F. If we assume that this operator has r distinct eigenvalues, then the 
spectrum can be written as 

cr(F) = {m,... ,a r }, ( 8 ) 

where in general ctj £ C. With each eigenvalue oti we associate its multiplicity denoted by n ? ;■ 
Obviously, £T =1 n i = where N refers to the dimension of the Hilbert space related to the 
analyzed system. 

The generator of evolution in (7) is still too general to consider the problem of observability. 
There might be many approaches to specify this generator and make the problem solvable. In this 
section, in order to get a concrete result, we assume that the operator F is unitary, i.e. F* = F _1 . 

Then the explicit form of the generator of evolution L can be obtained by using the idea of 
vectorization [15]. One can get 

L = (T _1 ) t <g>F + F T <g>F _1 - 2 Iat ® Hat. (9) 

For the generator of evolution presented in(9) we can formulate the following theorem and prove 
it. 

Theorem 2.1. The index of cyclicity for the generator of evolution given by (9) is expressed by 

r 

g = max{^2ni,5 1 ,...,5 p }, (10) 

i=l 


where 5k is defined as follows 

r—k 

4 = 2 y~] njUj+k (11) 

i —1 

and p = if r is odd or p = of r is even. 

Proof. One can instantly notice that the operator L is hermitian (self-adjoint), which implies that 
the algebraic multiplicities of its eigenvalues are the same as the corresponding geometric multiplic¬ 
ities. Therefore, we can focus only on determining the algebraic multiplicities. In order to do that 
one can realize that transposition does not change the spectrum of an operator. This observation 
leads to the conclusion that the eigenvalues of L are exactly the same as the eigenvalues of the 
operator L 7 which has the form 

L 7 = F <S> F 1 + (F ® F 1 ) 1 — 2 I 7 V ® I/v • (12) 

It can be noticed that the spectrum of F ® F is given by 

cr(F(g>F _1 ) = {mj € C\mj = —; i,j = l,...,r} (13) 

°Lj 

To each eigenvalue pij corresponds its multiplicity niUj. 

Bearing in mind equation (12), the spectrum of the generator L can be written as 

Q> • Q/- 

o-(L) = {A ij G M; Ay = — + — - 2; i,j = 1,... ,r} (14) 

at aj 

And multiplicity of A ij is equal mrij. Based on the equation (14) it is not possible to uniquely 
express the multiplicities of A ij by n t . But a specific case can be discussed when the eigenvalues 
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a±,... ,a r are successive terms of a geometric sequence. For clarity it can also be assumed that 
|cki| > |o: 2 1 > • • • > \ot r \- Let us denote 


where i = 2,... ,r-l. 

It is convenient to present the eigenvalues of L in a form of a matrix [A,j]: 


0 

9 + |-2 

9 2 + £ - 2 

9 + |-2 

0 

9 + |-2 

9 +^-2 

9 + |-2 

0 

9 3 + ^-2 
q 2 + X - 2 
9+|-2 

' ' ' 9 + 2 

••• q r ~ 2 + -^~2 

••• 9 r “ 3 + ^-2 


9 3 + ^-2 

9 2 + ^-2 

9 + |-2 

0 

n r ~ 4: _l_ 1 9 


f- 1 + - 2 

q r ~ 2 + ^=7 - 2 

q r ~ 3 + ^ - 2 

q r ~ 4 + ^ - 2 

0 

(16) 

It can easily be observed that the 

above matrix is 

symmetric and the eigenvalues that lie on 


the same diagonal are equal. Therefore, let us introduce denotation 5 k which will refer to the total 
multiplicity of all eigenvalues that lie on both k th diagonals (where k = 0 for the main diagonal). 
According to this definition the following equality holds 


r—k 

$k = 2 ^ ^ TljTli -|_fc (17) 

2=1 

It is sensible to consider only the cases for k £ {1, where p = if r is odd or p = if r is 
even, because for k > p the obvious inequality holds: 

r 

^2^ > 5 k . (18) 

2=1 

Therefore, the index of cyclicity for the analyzed generator can be expressed as 

r 

r] = max(^2nf,8i,...,5 p }, (19) 

2=1 

which ends the proof. □ 

In this section we have showed that under specific assumptions it is possible to calculate in a 
simple way the index of cyclicity for generators of evolution that consist of two Lindblad operators. 
Interestingly, the obtained final result is exactly the same as the index of cyclicity for the generator 
given by L p =-\[H, [H,p \] [9]. 

3 Superposition of double commutators as the generator of evo¬ 
lution 

Let us first formulate the problem which will be analyzed in this section. Now we assume to have an 

operator F £ B(H) which is hermitian, i.e. F* = F, and, as it was mentioned before, dimTL = N. 

Moreover we assume to have a set of matrix polynomials {/i,..., /at} ( fi : Mjy(C) Mjy(C)) and 
the degree of each matrix polynomial is not greater that N. 
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Then we shall define opeators F k in the way 

F k :=fk(F) for k=l,...,N. (20) 

Now we shall proceed to analyzing another model of evolution of open quantum systems. We 
propose to consider observability criteria for a generator of evolution given by 

N 

L P = J2 Lk ([■ F k p , F* k } + [F k , pF* k \) , (21) 

k =1 


where 7 ^ > 0 for k = 1 , ..., N. 

One can notice that we have F£ = F k because F* = F. Thus the evolution equation introduced 
in ( 21 ) can be expressed by means of double commutators 

N 

L p = — E 7 k[Fk,[ F k,p]\- (22) 

k= 1 

If one wanted to consider this problem in general, one would observe that it is impossible to 
reach a specific results due to a large number of parameters (one should bear in mind that each 
matrix polynomial f k can in general depend on N +1 coefficients). Therefore, in order to be able to 
draw concrete conclusions concerning conditions for quantum tomography, we propose to analyze 
the case when the matrix polynomials {/ 1 ,..., /jv} are defined in the following way 

f k (F) :=F k , for k = l,...,N. (23) 

This assumption allows to rewrite the evolution equation in the form 

N 

Lp = -^ k [F k ,[F k ,p}]. (24) 

k= 1 

Before introducing the main theorem we shall assume that the spectrum of the operator F has 
the following structure: r different eigenvalues denoted by {cci, • • • , a r } (oti € M) and with each a* 
we associate its multiplicity denoted by rij. As F in this model is assumed to be Hermitian, the 
geometric and algebraic multiplicities are equal. Therefore, there is no need to differentiate between 
them. 

For the generator of evolution given by (24) we shall prove the following theorem concerning 
the index of cyclicity. 

Theorem 3.1. The index of cyclicity for a quantum system with the generator of evolution given 
by (24) can be expressed by one of the three formulas below. 

1. such that 72^-1 / 0, then 


V = J2 n l ( 25 ) 

2=1 

If l 2 d-i = 0 for all d such that d € N and (2 d — 1 ) < N then we may have two situations: 

2. if r = 2u + 1 for u G N then 

2u-\-l u 

V = ^ ^ ^2 ^ ^ ^ TljTl2u-\-2—%•> 


( 26 ) 
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3. if r = 2 u for u € N then 

2u u 

3 = £ nf + 2 £ njn2u+i-i- (27) 

2=1 2=1 

Proof. To start the proof one needs to notice that the explicit matrix form of the operator L is 
given by 

N 

L = — £ 7fc ( F 2k ® Ijv - 2F fc ® + I N ® , (28) 

fc=i 

where the transposition that should appear according to vectorization equality has been skipped 
due to the fact that it does not influence the spectrum of L. 

One can then observe that the spectrum of the operator L is given by 

N 

(j(L) = {A ij <E M; A ij = - £ Tk( a i ~ for i,j = 1,.. ., r}. (29) 

k =i 

Important observation in this analysis is that A ij = 0 for a\ = a 3 . Therefore the total multiplicity 
of 0 as an eigenvalue of L (denoted by tuq) is at least equal to 

r 

-o = £n 2 . (30) 

2=1 

The question is whether this multiplicity can be greater under some additional assumptions 
concerning the structure of the spectrum of F. 

To prove point 1. of the theorem one needs to take into account the condition given in there, 
which means that there exists at least one non-zero 7 *, with k being odd, i.e. k = 2d — 1 for some 
d £ N. One can notice that with this condition \ t] = 0 only for a 3 = aj. Moreover in this case it 
is not possible to assume any symmetry between the eigenvalues of F that can lead to the situation 
in which the total multiplicity of some eigenvalue of L will be greated that the total multiplicity 
of zero. Therefore, the index of cyclicity in this case is the same as the (30), which completes the 
proof of the point 1 . 

In points 2. and 3. of the theorem 3 we consider the case when 7 *. for every k being odd is equal 
0. Under this assumption the generator of evolution can be simplified and written in the form 


L p 


If] 

£72 k [F 2 k ,[F 2 k ,p}}, 


(31) 


where symbol [x\ denotes the integer part of x. 

Consequently, the spectrum of the operato L can be written as 


o-(L) = {A i:j € M; Xij = - £ 72 k(ct 2k ~ a 2k ) 2 , for i,j = 1, - - • ,r}. (32) 

k =1 

Now we shall analyze two cases. 

1. Case for r = 2u + 1, where u £ N. 

To get a specific result we shall assume that the eigenvalues of F constitute an arithmetic 
sequence such that 

ot\ = uc where c £ R (33) 
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and 


cti+i — a.i = — c for i = 1, 


Bearing in mind this assumption one can notice that 


r — 1. 


(34) 


_ CpP o • — 0 

a 2w+2-i — 


(35) 


where i = 1, 2,... , u, u + 2, u + 3,..., 2u + 1 and k = 1,... , [^]. Equality (35) implies that 
\{ 2 u+ 2 -i) = 0, which means that the total multiplicity of zero as the eigenvalues of L can be 
expressed as 

2u+l u 

V = ^2 n * 2 + 2 X] n i n 2u+2-i , (36) 

2=1 2—1 


which is also the index of cyclicity of the generator of evolution and therefore the proof is completed. 
2. Case for r = 2 u, where u € N. 

To get a specific result we shall assume that the eigenvalues of F constitute an arithmetic 
sequence such that 


uc + (1 — i)c for i = 1,..., u; 
uc — ic for i = u + 1,..., 2 u. 


(37) 


One can easily notice that if the spectrum of F is defined in such a way, the following equality 


holds for k 6 {1,..., [4^]} 


a: 


2k 


— a. 


,2k 

222+1 — 2 


= 0 , 


(38) 


where i = 1,... ,2 u. This equailty implies that Aj( 2 u+i-i) = 0. Therefore the total multiplicity of 
zero as the eigenvalue of L is equal 


2u u 

V = n i + 2 n i n 2u+l-i , (39) 

2=1 2=1 

which completes the proof. □ 

In this section it has been proved that if one makes specific assumptions about the character of 
evolution of an N-level quantum system, one can obtain concrete results concerning conditions for 
observability. 

4 Conclusion 

In this paper we have formulated two evolution models of IV-level open quantum systems and 
we have applied to them the stroboscopic approach to quantum tomography. By making specific 
assumptions it was possible to determine for the analyzed generators of evolution the index of 
cyclicity, which physically refers to the minimal number of distinct observables that are necessary to 
perform quantum tomography. Because it focuses on determining the optimal criteria for quantum 
tomography, the stroboscopic approach has potential applications in experiments. This article 
contains two results connected with the criteria for quantum tomography, however the number 
of problems that can be solved with the stroboscopic approach is unlimited. One can treat the 
content of this article as a demonstration of possible applications of the stroboscopic tomography to 
selected evolution models. Consequently, the methods presented here can be used by one to solve 
more specific problems that one encounters in their research. 
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